Viewing the Proton Through "Color" -Filters 
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While the form factors and parton distributions provide separately the shape of the proton in 
coordinate and momentum spaces, a more powerful imaging of the proton structure can be obtained 
through phase-space distributions. Here we introduce the Wigner-type quark and gluon distribu- 
tions which depict a full-3D proton at every fixed light-cone momentum, like what is seen through 
momentum( "color" )-filters. After appropriate phase-space reductions, the Wigner distributions are 
related to the generalized parton distributions (GPD's) and transverse- momentum dependent parton 
distributions, which are measurable in high-energy experiments. The new interpretation of GPD's 
provides a classical way to visualize the orbital motion of the quarks, which is known to be the key 
to the spin and magnetic moment of the proton. 



What is the shape of the proton? The textbook an- 
swer is that it is round 1]. Indeed, the proton is a spin- 
1/2 particle, and the strong interactions are invariant 
under time-reversal symmetry. Hence there cannot be 
any dipole, quadrupole or higher-multiple deformations, 
and the quark and gluon density distributions in space 
have to be spherically symmetric. It is well known that 
the Fourier transformation of the proton's electric form 
factor yields the spherical charge distribution when the 
relativistic recoil effects are neglected PJ . 

The classical answer is correct, but might not be com- 
plete. The quarks and gluons are not static inside the 
proton, they are moving relativistically to produce the 
magnetic moment and part of the angular momentum 
(spin) of the proton. The orbital motion has an axial, 
rather than spherical, symmetry when the proton spin is 
polarized along a fixed direction (3j. To be sure, the Feyn- 
man distributions do provide a momentum space picture 
of partons, but there is no spatial density information 
involved — they are just correlation functions in the posi- 
tion space. To gain a more complete understanding of the 
proton's structure, it is useful to know the quarks' and 
gluons' joint-position-and-momentum (or phase-space) 
distributions. Indeed, in classical physics, a state of a 
particle is specified by knowing both its position x and 
momentum p. In a gas of classical identical particles, the 
single-particle properties are described by a phase-space 
distribution f{r,p) representing the density of particles 
at phase-space positions (r , p) . Time-dependence of the 
distribution is governed by the Boltzmann equation. 

In quantum mechanics, however, the notion of a phase- 
space distribution directly contradicts one of its funda- 
mental principles: the momentum and position of a par- 
ticle cannot be determined simultaneously. Nonetheless, 
over the years, physicists have introduced various quan- 
tum phase-space distributions which reduce to f{r,p) in 
the classical limit. These distributions have proven ex- 
tremely useful in wide-spread areas like heavy-ion col- 
lisions, quantum molecular dynamics, signal analysis, 
quantum information, non-linear dynamics, optics, im- 
age processing, etc. One of the most frequently used is 



the Wigner distribution |2| 

W{x,p)^ J i;*{x-r]/2)i;{x + 7j/2)e'P''dri , (1) 

where ipi^) is a quantum wave function (h = 1). When 
integrating out the coordinate x, one gets the momentum 
density |'!/'(p)P, and when integrating out p, the coordi- 
nate space density |'i/;(a;)p follows. For arbitrary p and x, 
the Wigner distribution is not positive definite and does 
not have, strictly speaking, a probability interpretation. 
However, there is no doubt that it contains much more 
information than the usual quantum mechanical observ- 
ables. Modifed Wigner distributions with positivity are 
known in the literature and are better interpreted as den- 
sities, but we leave them for a further study. 

In this paper, we introduce the quark and gluon 
Wigner-type phase-space distributions in hadrons such as 
the proton. After certain phase-space reductions, these 
distributions reduce to the generalized parton distribu- 
tions (GPD's) and the transverse-momentum dependent 
parton distributions which are measurable in high-energy 
experiments. [For reviews of these distributions and ex- 
periments to measure them, see 0, 0|-] Their relation 
to the GPD's reveals a new interpretation of the latter — 
the Fourier transformation of GPD's yields full-3D quark 
and gluons images of the proton at every fixed Feynman 
momentum x, like what's seen through momentum (or 
better yet "color" or x-) filters. The skewness parameter 
^ in GPD's now has a natural interpretation: a measure 
of the proton deformation along the direction selected by 
high-energy probes. From the sign variation of the im- 
ages, we might have a practical way to determine how 
"classical" the quarks and gluons are in the proton. 

Disregarding the renormalization scale fi, the GPD's 
depend on three kinematic variables x (Feynman mo- 
mentum), ^ (skewness), and (t-channel momentum 
transfer, also called t). Introduced as quantum interfer- 
ence amplitudes, their original interpretation as parton 
and their angular momentum densities was limited to 
t = ^ = 0|2l, Izl- In a seminal paper by M. Burkardt 
the probability interpretation was successfully extended 
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to the GPD's at t 0: the Fourier transformation of the 
parton densities in the two-dimensional transverse plane 
(impact-parameter space). A number of subsequent pa- 
pers have studied the physical significance of the ^ ^ 
case , without seeking a classical or quasi-classical in- 
terpretation. The present work expands the 2D-pictures 
of Burkardt along the longitudinal direction, relating the 
GPD's to full-3D distributions at a fixed light-cone mo- 
mentum, interpretable as densities in the classical limit. 

We begin by introducing the Wigner distribution op- 
erator for quarks in QCD 



Wr{r,k)= / ^{f-i]/2)T^{r + i]/2)i 



(2) 



where r is the quark phase-space position and k the 
phase-space four-momentum. P is a Dirac matrix defin- 
ing the types of quark denstities. VE" is a gauge-invariant 
quark field in non-singular gauges (gauge potentials van- 
ish at the spacetime infinity [lOj), 

*(?7) = exp (^-ig J dXn ■ A{Xn + 77)^ V('?) , (3) 

where n'^ is a constant four- vector. Thus VVr is gauge- 
invariant but depends on a choice of n^^. The actual 
Wigner distribution can be define as the expectation 
value of Wr in the hadron states. For example. 



Wr{r,k) = 



1 

2M 
1 

2M 



(27r)3 
d^q 
(27r)3 



q/2 W(f, k) - q/2 



(4) 



^^■Uq'/2 WiO,k) ~q/2 



where M is the hadron mass and the states are normal- 
ized relativistically. The initial and final hadrons are 
taken with different center-of-mass momenta, otherwise 
translational invariance results in a trivial dependence 
on the phase-space position r For non-relativistic 
quantum systems, the above definition readily reduces to 
the known Wigner distribution. For relativistic systems 
such as the proton, the above interpretation as a density 
neglects relativistic effects when the momenta are com- 
parable to the proton mass Fortunately, since the 
proton mass is of order 1 GeV and charge radius 0.86 fm, 
the relativistic correction is about %10. 

For the proton, because of the high energy involved in 
experiments and because of the color confinement, the 
k~ = (fc" — k^)/^/2 energy distribution is difficult to 
measure, where the z-axis refers to the momentum di- 
rection of a probe. Moreover, the leading observables 
are associated with the "good" components of the quark 
(gluon) fields in the sense of light-cone quantization , 
which can be selected by P = 7+, 7^75, or a^-^ where 
7^ = (7°-|-7^)/-\/2. The direction of the gauge link, n'', is 
then determined by the trajectories of high-energy par- 
tons traveling along the light-cone (1,0,0,-1). There- 
fore, from now on, we restrict ourselves to the reduced 



Wigner distributions integrating over k 



Wr{r,k)^ 



dk- 

W) 



rWr(r,fc) 



(5) 



with a light-cone gauge link. Eq. ||SJ defines the most 
general (master) phase-space quark distributions rele- 
vant in high-energy processes. Unfortunately, there is no 
known experiment at the moment capable of measuring 
them in the full 6-dimensional phase space. 

However, certain phase-space reductions of the distri- 
butions are measurable through known processes. For 
example, integrating out the transverse momentum of 
partons, we find 



/^(r,fc^ 



J^Wr (r,k) 
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2M 
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(27r)3 J 27r 

X {q/2 |*(-rr/2)P*(7r/2)| 



q/2) (6) 



The matrix element under the integrals is what defines 
the GPD's. More precisely, if one replaces fc+ by Feyn- 
man variable xP+ (P+ = Eq/^/2, proton energy Eg — 
^ Af 2 -I- (7^/4) and r/^ by A/P+, the reduced Wigner dis- 
tribution becomes the Fourier transformation of the GPD 



frir,x) 



1 

2M 



d^q 
(27r)3 



^«"Fr(a;,e,t) . (7) 



In the present context, the relation between kinematic 
variables are ^ = qz/{2Eq) and t = —q^- Taking P 
V27^, the corresponding GPD has the expansion 



^e'^^ (q/2 V'(-An/2)/:V27+V(An/2) 
27r \ 

H{x, t)U{q/2)V2j+U{-q/2) 
+ E{x,tt)Um)'-^U{-q/2) , 



q/2 



(8) 



where C is the shorthand for the light-cone gauge link. 

The distribution /^+ (r , x) can be interpreted as the 3D 
density in the rest frame of the proton, for the quarks 
with a selected light-cone momentum x. Integrating over 
the z coordinate, the GPD's are set to ^ ^ = 0, and 
the resulting two-dimensional density f^+{fj_,x) is just 
the impact-parameter space distribution in the Burkardt 
paper ^. 

To understand better the physical content of the den- 
sity, let us examine its spin structure. Working out the 
matrix element in Eq. (|SJ), 



1 

2M 



P^+(x,C,i) = [H{x,^,t)-TE{x,^,t)] 



+ iisxq)'—[H{x,C,t)+Eix,^,t)] 



(9) 
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where r = (1 + ^)g^/4M^ The first term is indepen- 
dent of the proton spin, and the related quark density is 



p+{r,x) 



(27r)3 



e-'«-^Xa;,e,t)-r£;(a;,e,0] • (10) 



Just hke charge density, it is not positive-definite at any 
X. The second term depends on the proton spin, the 
corresponding density is physically the 3rd-component 
of a current 



^ [Hix,C,t)+Eix,^,t)] 



2M 



(11) 



The i7-term current generates the proton's Dirac mo- 
ment, and the iJ-term generates a convection current 
due to the orbital angular momentum of massless quarks 
and vanishes when all quarks are in the s-orbit[l3. The 
physics in separating /+ into /?+ and can be seen from 
the Dirac matrix 7+ selected by the high-energy probes, 
which is a combination of time and space components. 
Because the current distribution has no spherical sym- 
metry, the quark density seen in the infinite momentum 
frame, /9+ -I- j^, is deformed in the impact parameter 
space 13]. This is the kinematic effect of Lorentz trans- 
formations. 

We can also study the a;-moments of the quark and cur- 
rent densities. If we integrate over x directly, we obtain 
the electric charge and current densities, 



P+i^ ^ I dxp+{r,x) 



(27r)3 
d^q 



1 



{2n) 



^e-^^-\{sxq) — {F,{q^)+F^{q^)) 



Notice that p'^j^ is not spherically symmetric due to the ^ 
dependence in r, reminiscent of the 7"'"-probe. The elec- 
tric current density has a donut shape and is responsible 
for the magnetic moment of the proton Q . Indeed, if we 
calculate the proton's magnetic moment using the clas- 
sical formula, p = (1/2) / r x j, the right answer follows. 
The current generated by the ^2(9^) term is directly re- 
lated to the orbital motion of the quarks Q| . 

Forming the second moment of a;, we have the mass 
distribution of the proton from the quarks , 



dxxp^ if, x) 
(27r)3 



-iq 



- [{A,{q')+eCq{q')) 



'{B,{q^)-eCq{q^)] 



(12) 



where Aq{q^), Bq{q'^) and Cq{q^) are the quark part of 
the proton's energy-momentum form factors 



larkpa 



A{q^)P''^-^ + B{q'^'^ '^"^ 



2M 



+Ciq'){q'^q-'-g^''q')]UiP) . (13) 

Here the energy-momentum tensor T**" is symmetric and 
conserved, and contain both quark and gluon contribu- 
tions. The normalization of the form factors ^(0) and 
B{0) are well-known: A{0) = 1 is related to the energy- 
momentum conservation and B{0) = is a result of the 
angular momentum conservation Il7l| . 

The energy-momentum form factors can in principle 
be measured in elastic graviton-nucleon scattering. Here 
we use this fact to motivate the physical interpretation of 
p™(r) as the proton's mass distribution. Gravitions have 
spin 2, and an on-shcU massless graviton has helicity A = 
±2. An off-shell (massive) graviton can have additional 
helicity states A = ±1 and 0. Consider the graviton- 
proton scattering in the Breit frame. Assume the helicity 
of the initial proton is 1/2, then the final proton has 
helicity 1/2 — A. Since the proton has spin- 1/2, only off- 
shell gravitons of helicity 1 and contribute. A graviton 
with helicity 1 and acts very much like a photon, and 
we can mimic the Sachs's electromagnetic form factors 
by defining the gravitational "electric" (or mass) form 
factor 



G^{q') = A{q^) + 



4Af 



(14) 



proportional to the helicity non-fiip graviton-nucleon 
scattering amplitude [l^ UJj . If one ignores the defor- 
mation due to the probes (^ = 0), p™(r) in Eq. H12|l is 
just the Fourier transformation of the above form factor. 

The x-moment of the quark current is the momentum 
density in the proton 



M 

— I dxxj^ 
d^q 



dxxj'^ (r , x) 
^e~'^-^i{sx q)' 



(2^) 

^^[Me) + Bq{q^)\ . (15) 

The combination in the integrand is the gravitational 
"magnetic" (or angular momentum) form factor. 



G^A<l')^A{q') + B{q') 



(16) 



proportional to the helicity-flip graviton-proton scatter- 
ing amplitude 0, 0| . Knowing the momentum density 
as a function of r, it is simple to calculate the quark con- 
tribution to the angular momentum (spin) of the proton 
according to / f x . The answer is {Aq{Q) + Bq{0))/2 
A contribution from the gluons can be obtained in a 
similar way. 
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Let us remark that gravity plays a dominant role at 
two extreme distance scales: the cosmic scale and the 
Plank scale. At the atomic and subatomic levels, with the 
exception of the motion of these particles under strong 
gravitational fields, gravitation has little effects because 
of its extremely weak coupling. Thus while mass distri- 
bution and moment-of-inertia are important concepts in 
classical mechanics, they have rarely been discussed for 
microscopic quantum systems such as the proton. Thus, 
it is very interesting that the GPD's provide gravita- 
tional form factors without actual gravitational scatter- 
ing 1^, 0]. A study of the gravitational Sachs's form 
factors (and gravitational radius) in chiral perturbation 
theory can be found in A related coordinate-space 
interpretation of the gravitational form factors, in par- 
ticular, the stress density T*-' , can be found in Ref. [20j. 

If integrating over r in the Wigner distributions in 
Eq. jSJ, one obtains the transverse- momentum depen- 
dent parton distributions. There is a lot of interest- 
ing physics associated with these distributions which has 
been discussed recently. For instance, in a transversely 
polarized proton, the quark momentum distribution has 
an azimuthal angular dependence |^ Hol |2]| . The so- 
called Siver's function can produce a novel single-spin 
asymmetry in deep-inelastic scattering. We will not pur- 
sue this topic here, except emphasizing that they have 
the same generating functions as the GPD's. 

The phase-space distributions in the transverse coordi- 
nates have also been used in small-x physics 22], where 
the parton physics becomes truly classical at high den- 
sity 

We end the paper by making a few remarks. First, as 
we have mentioned before, the Wigner distributions are 
quantum phase-space distributions which contain much 
more information of a quantum system than conventional 
observables. Originally introduced as a theoretical con- 
struction, it is surprising that these distributions for the 
proton are directly measurable through high-energy pro- 
cesses. Second, the Wigner distributions are more than 
just quantum interference amplitudes; they have a clear 
physical interpretation in the classical limit. This endows 
the GPD's a new physics feature which has not been rec- 
ognized before. QCD in the limit of large number of 
colors might be a case in which the Wigner distributions 
become truly classical densities. Finally, the extent to 
which the Wigner distributions in the proton can be re- 
garded as classical densities may be judged once the phe- 
nomenological GPD's is available from experiment data. 

I first heard a possible connection between GPD's and 
Wigner distribution from A. Bclitsky, and this work is 
inspired by him and T. Cohen. I also thank M. Burkardt 
for many helpful discussions about his impact-parameter 



space interpretation of the GPD's, and M. Diehl for crit- 
ical comments. This work was supported by the U. S. 
Department of Energy via grant DE-FG02-93ER-40762. 
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